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Abstract 

We study the extended hydrodynamics proposed by PhiUp Rosenau [Phys. Rev. A 40, 7193 
(1989)] in the context of a regularization of the Chapman-Enskog expansion. We are able to prove 
that shock waves appear in finite time in Rosenau's extended Burgers' equation, and we discuss the 
physical implications of this fact and its connection with a possible extension of hydrodynamics to 
the short wavelength domain. 
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The Boltzmann equation is one of the most fundamental equations in Nonequihbrium 
Statistical Mechanics. This equation describes the dynamics of a rarefied gas, taking into 
account two basic processes: the free flight of the particles and their collisions. Due to the 
difficulties that a direct treatment of this equation implies, a reduced description of Boltz- 
mann equation is one of the major problems in kinetic theory. Equations of hydrodynamics 
constitute a closed set of equations for the three hydrodynamic fields: local density, local ve- 
locity, and local temperature. These equations can be derived from Boltzmann equation by 
performing the Chapman- Enskog expansion Q . This expansion is a power series expansion 
in the Knudsen number, that is the ratio of the free mean path between the macroscopic 
length. The first order of the expansion yields Euler equations, while the second order yields 
Navier-Stokes equations which in the case of an incompressible fiuid read: 

dtv + (v • V)v = -Vp + /iV^v 

V-v = 0, (1) 

where /i represents the viscosity of the fiuid, and is of the order of the Knudsen number. The 
next order in the Chapman- Enskog expansion yields the Burnett equations of hydrodynamics 
which are, unfortunately, invalid. To see this more clearly consider the viscous part of the 
Chapman- Enskog expansion: 

e(^oVV + eViW + ...), (2) 

where n = efio, and e is the Knudsen number. The Burnett order implies the presence of the 
biharmonic term proportional to V^, that causes an unphysical increase in the number of 
boundary conditions and rends the equilibrium unstable, among other undesirable effects. 
While the Navier-Stokes equations give very accurate results in many domains, they usually 
fail when applied to predict the short wavelength properties of the fluid, like, for instance, 
the propagation of ultrasounds within the fluid. This makes very useful to develop a higher 
order description of the fluid, while the Burnett order has proven itself less accurate than the 
Navier-Stokes order. This problem was partially solved by Philip Rosenau in his influencing 
article of 1989 P]. The idea was to regularize the Chapman- Enskog expansion using a very 
original comparison. First consider the power series expansion: 

-^ = l + z^ + z^ + ..., (3) 
2 



where 2; is a complex number which modulus fullfills \z\ < 1. Assuming that e is small enough 
and taking into account the power series Q suggest that we can recast expansion ^ into 
the form 



1 — e^m^V^ ' 

where rn? = — , and this operator is to be interpreted in the Fourier transform sense: 



(4) 



Mo ' 



1 — e^m^V^ / 1 + e^m^/c^ 
This idea was originally proposed in the__context of random walk theory jsj, and has been 



used within this context in latter works 

While this regularization of the Chapman-Enskog expansion seems to be a proper exten- 
sion of hydrodynamics in the linear regime Q , its effect on the full nonlinear hydrodynamics 
is not so clear. This is due to the analytical difficulties that a mathematical treatment of 
the Navier-Stokes equations implie. However, it is useful to study some toy models to win a 
deeper understanding of hydrodynamics; to this end was developed a one- dimensional model 
for hydrodynamics: the Burgers' equation 

dtu + udxU = fxdlu. (6) 

In the same spirit, Rosenau considered the regularized Burgers' equation, arguing that an 
understanding of this model would clarify the effect of the regularization of the Chapman- 
Enskog expansion on the nonlinear hydrodynamics. The rest of this work is devoted to prove 
the appearence of shock waves in finite time in Rosenau's regularized Burgers' equation, and 
to analize the physical implications of this fact. 
Rosenau's extended Burgers' equation reads: 

dtU + ud^u = fi _ ""^^^ M, (7) 

where we have set, without lost of generality, fii/ fiQ = 1. To prove shock wave formation we 
will exploit the analogy between viscous Burgers' equation (the inviscid Burgers' equation 
is obtained just by setting /i = 0) and the Keller-Segel system |5]: 

dtv = fid^v + d^ivd^w), (8) 
dlw = -V. (9) 
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Note that we recover viscous Burgers' equation performing the substitution u 
system fj8l9|) . Consider now the following modified Keller-Segel system: 



dtv 



1 - eW, 

V. 



V + d^ivd^Wj, 



dxW in 

(10) 
(11) 



We can recover Rosenau's extended Burgers' equation by performing again the substitution 
u = dxW in this last system. We will consider homogeneous Dirichlet boundary conditions: 
v\dn = w\q^ = 0, where Vt is the closed interval VL = [— L, L]. From system (jlUllip we get: 



1 d 
2di^ 



f 92 

VVtdx = /i / V ^ 



vdx 



X 

vdxwdxvdx + / v^dx (12) 
n Jn 



Now, we are going to estimate all the terms appearing in the right hand side of this equation. 
Integrating by parts the second term in the right hand side of Eq. lfT^ : 



vdxwdxvdx = v'^dxw\Q^ — / dxvdxwvdx 



vd^wvdx, 



(13) 



that implies 



vdxwdxvdx = —7: V d^wdx 



- / v^dx 
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The first term in the right hand side of Eq. ()12|) can be estimated as follows: 

5! 



(14) 



vdx < 
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1 - e^dl 



i2 



< 



dx < \ \v\ 



L2(n) 



1 - e'^dl 



(15) 



where we have used Holder's inequality (see below). By performing the shift of variables 
y = x/e, we get: 
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1 - e292 



L2(0) 



.(3/2) 



1-^2 

y 
N 



< 



L2(n/e) 



.(3/2) lrML2(C/e) ' 



(16) 



where N = [5^(1 — Let us clarify a bit this last step. We have used the fact that the 

operator V^(l — V^)~^ is bounded on every space, with I < p < oo. This means that 
we can assure that ||V^(1 — ^'^)^^ f\\Lp{n) — ^ \ \f\\Lp{n) ^'^^ every / belonging to Lp{Q) and 
a constant that does not depend on / (and thus is called the norm of the operator). 
This fact can be easily seen once one realizes that the Fourier transform of the operator 
V^fl — V^)~^ is a bounded function of the wavevector, and a rigorous proof can be found 

n 

in p]. We can again shift variables x = ey to get: 

.<^.<^IMIL.. (IT) 



Finally, we can conclude our estimate as follows: 
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-vdx > - 



n 



1 - e^dl 



N M m2 . ^ 

—^mym- (18) 



Now we are going to estimate the third term in Eq. ()12p : 



v^dx = \\v\\l;,^... (19) 



L3(n) 

Holder's inequality reads (for a rigorous proof of Holder's inequality see Q): 

\fg\dx < \\f\\Lv(^n) ll^llL^n)' 

l<p,g<oo, 1 + 1 = 1. (20) 



Choosing = 1 we get: 



where C = . With this estimate we can claim that: 



ljf\dx<C\\f\\,,^^^, (21) 



1^11^2(0) = / v^dx<C\y\\^^^^^ 



\ (i/p) 



/ r \(2/3) 

C[J^^'dxj =C\\v\\l^^^, (22) 
where we have chosen p = 3/2 (and correspondingly g = 3). This implies that: 

\\v\\LHn)>D\\v\\L2^n)^ (23) 



where D = Therefore, we have the final estimate: 

j^Mll^^n) > ^ {M\%(n)y ' ^ - ^ , (24) 

where A, i? > are constants, A = and B = We are thus going to study the 
dynamical system: 

^ = - Bx. (25) 

dt ^ ' 

This system has two fixed points, x = and x = (B/A)^ > 0. A linear stability analysis 

reveals that the positive fixed point is linearly unstable, meaning that every initial condition 

Xq > {B/AY will stay above this value for all times. Further, we know that the solution will 

grow without bound in this case, so we can claim the existence of two constants, < oo 

and < Co < v4, such that Ax'^l'^it) - Bx{t) > Cox^/'^(t) for every t > to. This implies that 

^ II i|2 (w i|2 V'^'^ 

for t > to, and for an adecuate initial condition. Solving this equation gives: 

H-MUn) > I == (27) 

for t > ti > to, and for an adecuate initial condition. And every adecuate initial condition 
must fuUfiU 

\H-M\lHn)>^^\^\ 
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+ ^ll^(-,0)||^.(^)+ l^||i;(-,0)||i.(^), (28) 

like, for instance, v{x,0) = (a;^ + 6)~^^^ — (L^ + 6)~^^^ and 6 small enough. Thus we are 
finally led to conclude that the system does blow up in finite time. If we recover v = —d^u 
we see that the first spatial derivative of u becomes singular in finite time. This means that 
the solution to equation Eq.Q develops a shock wave in finite time (or what is the same, a 
discontinuity in the fiow appears), in contrast to the viscous Burgers' equation and analo- 
gously to the inviscid Burgers' equation dtv = —vdxV ^ . The inviscid Burguers' equation is 
a one- dimensional model for the Euler equations, while the viscous Burguers' equation sim- 
ulates the Navier-Stokes equations. This suggests that the regularizating procedure implies 
a return to a lower order in the Chapman-Enskog expansion. 
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It was already argued by Rosenau that this kind of regularization of the Chapman- Enskog 
expansion was only valid in the linear regime, while nonlinear terms might be present in 
the full nonlinear hydrodynamics. These terms are expected to have a deep impact on 
the dynamics of the fluid, the reason being as follows. Whether or not the Navier-Stokes 
equations become singular in finite time is still unkown and it is actually one of the most 
important open problems in Mathematics. What one would expect from a physical point of 
view is that these possible divergences smooth out if we look closer to the fluid taking into 
account higher order terms in a (complete) regularized Chapman-Enskog expansion. What 
we have seen in this work is that the linear regularized theory is able to convert a smooth 
solution into a singular one, so one would expect that a regularized Navier-Stokes equation 
of the form: 

^t^r + (v ■ V)v = -Vp + fi- 5^^v 

V • V = (29) 

is less regular than the original Navier-Stokes equation. We expect the presence of this 
nonlinear terms to regularize enough this equation that one would be able to prove global 
existence in time of the solution, and this way give a precise physical meaning to the possible 
divergences arising in the original Navier-Stokes equation. 
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